Anomalous slow fidelity decay for symmetry breaking perturbations 
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Symmetries as well as other special conditions can cause anomalous slowing down of fidelity 
decay. These situations will be characterized, and a family of random matrix models to emulate 
them generically presented. An analytic solution based on exponentiated linear response will be 
given. For one representative case the exact solution is obtained from a supersymmetric calculation. 
The results agree well with dynamical calculations for a kicked top. 
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Sensitivity to perturbations as measured by fidelity 
decay has received a great deal of attention both in 
the context of quantum and classical dynamical sys- 
tems 0, 12- S Q and as a benchmark for the stability 
of possible quantum information devices dEl. Fidelity 
may be described as the cross-correlation function be- 
tween unperturbed and perturbed evolution of a quan- 
tum or classical wave system. Linear response theory has 
been particularly successful in describing fidelity decay, 
by relating it to the correlation decay of the perturbing 
operator in the interaction picture 0, 0, El • For chaotic 
systems a recently introduced random matrix model 
is in excellent agreement with experiments [H . 

It has been shown, that so-called residual perturba- 
tions that have vanishing diagonal matrix elements in 
the eigenbasis of the unperturbed Hamiltonian, lead to 
very slow fidelity decay known also as quantum freeze 
of fidelity 0, [Tj- Fidelity freeze is a pure wave phe- 
nomenon without classical analogue In the field of 
quantum computation, it is known as dynamical decou- 
pling |13L Il4| . One can identify four different physical 
situations in which quantum freeze or anomalous slow fi- 
delity decay occurs. The case when the perturbation can 
be written as a time derivative (or commutator with the 
unperturbed Hamiltonian) is treated in Refs. 0, 0] . 

In this letter we treat the three other and indeed phys- 
ically most important cases: The first one corresponds 
to a perturbation which breaks an antiunitary symmetry 
(e.g. time-reversal) in an optimal way, meaning that the 
perturbation anticommutes with the antiunitary symme- 
try, e.g. switching on the magnetic field. The other two 
cases correspond to a mean field approach in which the 
diagonal part of the perturbation is moved to the un- 
perturbed Hamiltonian. We will consider unperturbed 
Hamiltonians, with and without an antiunitary symme- 
try. To obtain a generic understanding, all cases are 
considered in the framework of random matrix theory 
(RMT). We present general theoretical results in the lin- 



ear respons regime, as well as an exact analytical result 
obtained by a supersymmetry method for the antiu- 
nitary symmetry-breaking case. Our results display ex- 
cellent agreement with numerical experiments for quan- 
tum kicked tops. 

To study fidelity decay, we consider the perturbed 
Hamiltonian H = H + XV. If U(t) and U Q (t) are the 
unitary propagators under H and Ho, respectively, we 
define the fidelity amplitude as 

f(t) = (*o(OI*(*)> = (*(0)|E*>(-t) U(t)\9(0)) , (!) 
where #(i) = U(t) *(0) and # (*) = UMU!(0). Fidelity 



fa, 



is defined as F(t) = \f(t)\ 2 . In Refs. |2j-L|. it has been 
shown that within second order time-dependent pertur- 
bation theory (Born series) the fidelity amplitude, aver- 
aged over random initial states, can be expressed in terms 
of the two-point time correlation integral C(t): 



(f(t)) E = 1 - 4tt^ C(t) + 0(A 4 ) , 



(2) 



where (...)e denotes the average over random initial 
states, which are concentrated on a small energy interval 
that contains many levels Ne- For the fidelity ampli- 
tude, this averaging amounts to taking a restricted trace 
of the echo operator, in the eigenbasis of the unperturbed 
Hamiltonian Hq. Therefore, we may write: 



C(t) 



dt' / dt"(V(t')V(t")} 1 



(3) 



(V(t')V(t")) E = — 



E E Wo, 



2 2m.(E a —Ep)(t'—t" 



where the E a are the eigenvalues of H , and V(t) — 
Uo(ty V Uo(t) is the perturbation in the interaction pic- 
ture. The matrix elements V a p are taken in the eigenbasis 
of Hq. The primed sum runs over N E eigenstates of Hq. 
Often, the exponentiated version of Eq. |(SJ|, 



(f(t))E = e 



o(\ 4 



2 \2 



(4) 
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is able to describe the fidelity decay from t he p erturba- 
tive to the golden rule regime well 0, @, 0, EI ■ This is 
true, in particular, if the system is chaotic, such that the 
perturbation can be described by one of the random 
Gaussian ensembles (RMT approach) [l5j . 

We now consider the situation of a residual perturba- 
tion, i.e., one with vanishing diagonal, V aa = within 
the RMT models. Thus, we assume that the non-zero 
matrix elements of this perturbation are independent 
normalized Gaussian random variables with variance 



This yields a logarithmically slow decay of fidelity 



{\V aP \ 2 ) = 1-S afi 



(•5) 



where (...) denotes an ensemble average. For the pertur- 
bation matrix V, we consider three different ensembles, 
which fulfill Eq. (0): (i) an ensemble of imaginary anti- 
symmetric matrices, (ii) an ensemble of real symmetric 
matrices with deleted diagonal, and (iii) an ensemble of 
Hermitean matrices with deleted diagonal. 

The average of the correlation integral C(t) over any of 
those ensembles gives: 



<£(<)} = 5 



At' / dt"b(t") 



(6) 



where b(t) is the two-point spectral form factor of Ho. 
Here, as well as throughout the rest of the letter, we use 
units, where the Heisenberg time is equal to one. If Ne 
is sufficiently large, b{t) tends to a well defined smooth 
function (self-averaging); else we average over a random 
matrix ensemble for Hq as well. For Gaussian orthogonal 
(GOE) and unitary (GUE) ensembles, the form factors 
are given in Ref . [15| • Note that the term proportional to 
t 2 is missing as compared to the linear response result for 
a generic perturbation |2|. For what follows, we assume 
that the average (. . .) includes such a procedure. 

Let us first concentrate on the special case of unper- 
turbed Hamiltonians Hq taken from the GUE. Then, we 
have b{t) = max{l — t, 0}, and 



(C(t)> = C GUE (i) 



t < 1, 
t > 1 



(7) 



As a result, for times longer than the Heisenberg time, 
Eq. @ predicts the fidelity to "freeze" on a minimal value 



/plateau — 1 ~ 
6 



(8) 



Since the next correction term grows quadratically in 
time, (f(t)) = / p iatcau + C(A 4 t 2 ), we find that the plateau 
ends at a time of order t* ~ 1/A. 

Second, we choose Hq from the GOE. Also in this case, 
the integrals in equation © can be performed analyti- 
cally (see Ref. p|). Here we just give the the leading 
asymptotics for t I: 



CgoeW 



ln(2t) 



12 



0{t~ l Int) 



(9) 



(fit)) 



l-±[2 + ln(2t)] + 0(\H 2 ) 



(10) 



Both results for the plateau value of the fidelity ampli- 
tude, Eq. © and Eq. lO, follow from Eq. ©. They 
are thus valid for any of the three ensembles used for the 
perturbation, (i), (ii), and (iii). Indeed, a similar result 
could be obtained for the Gaussian symplectic ensemble. 

In the averages (F(t)) and |(/(i)}| 2 were shown to 
differ only in a term proportional to £ 2 , which is absent 
in the quantum freeze case. Thus, the eqs. © and © 
yield the linear response approximation for fidelity decay. 
Numerics indicate that (F(t)) « \(f{t))\ 2 , in general, in 
accordance with an argument given in . 

For long times and small perturbations, the fidelity am- 
plitude can be expressed in terms of level shifts. Within 
the second order stationary perturbation theory, its av- 
erage is then given by the Fourier transform of the level 
curvature distribution [l^j ] which was obtained analyti- 
cally in [lflj . The final result is surprisingly simple: 



(/(*)> 



tK^t) 
(1 + t) e- T 
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GUE 



it 



(11) 



where K\ is the modified Bessel function of first or- 
der. The GOE branch is valid for an unperturbed GOE 
Hamiltonian, and a perturbation matrix of type (i) or (ii). 
The GUE branch is valid for an unperturbed GUE Hamil- 
tonian, and a perturbation matrix of type (iii) . More de- 
tails on the asymptotic behavior of fidelity decay in situ- 
ations of quantum freeze will be published elsewhere . 

One should stress that diagonal elements of the pertur- 
bation vanish also in the presence of a discrete or contin- 
uous unitary symmetry R, of Hq, which anti-commutes 
with V, RV — —VR. However, it turns out that its effect 
on fidelity enhancement is less drastic than the predic- 
tions of Eqs. and JJjJ, because of the lack of corre- 
lations between different subspectra of Hq. As a result, 
the asymptotic growth of the correlation integral is linear 
C(t) oc t, for times before and after the Heisenberg time. 

For the case of Hq taken from the GOE and a purely 
imaginary antisymmetric perturbation the average fi- 
delity can be obtained exactly by supersymmetry tech- 
niques in the limit of large dimension N. The calcula- 
tion is technically much more involved than for the case 
of a GOE perturbation and will be presented else- 
where 18]. The result again is a VWZ-likc integral (see 
Ref. HaTEq. (8.10)) and is given by 



(f(t))=2 / du 



' dv 



V> 2 -v 2 ][(u + l) 2 -v 2 } 
[l + e{t 2 ~v 2 )] 



Max(CU-l) 

(t-u)(l-t + u) 
(v 2 - t 2 ) 2 

X [t(2u + 1 - t) + v 2] e -m^+l-t)-v*] 



(12) 
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FIG. 1: The complement of the average fidelity amplitude, for 
a weak perturbation, e = 5 x 10~ 5 . Part (a) shows the GOE 
case with a purely antisymmetric random Gaussian pertur- 
bation. Part (b) shows the GUE case, with an independent 
GUE perturbation with deleted diagonal. The exponentiated 
linear response approximations are plotted with short dashed 
lines (for comparison, the exact results for full GOE and GUE 
perturbations are shown with chain curves) . The results from 
time-independent perturbation theory are plotted with long 
dashed lines. In the GOE case, the exact analytical result, 
Eq. 112H is plotted with a solid line. The random matrix sim- 
ulations are plotted with points. 



The only difference to Ref. [lfj , where a GOE perturba- 
tion was considered, is the additional factor [l+e(t 2 — v 2 )] 
in the integrand, and a minus sign with the v term in 
the exponent, where in the GOE case there is a plus sign. 

Figure ^ shows the fidelity decay both for Ho taken 
from the GOE (a) and the GUE (b), for a small pertur- 
bation e = 5 x 10~ 5 . It compares different theoretical 
approaches. To this end, one minus the average fidelity 
amplitude is plotted on a double-log scale. The figure 
shows the exact result calculated from Eq. 1)12)) (in the 
GOE case only), together with the results from the ex- 
ponentiated linear response approximation J3j and the 
asymptotic expression 1)110. For comparison the fidelity 
decay for the case of a GOE perturbation 0] is shown 
as well. We see that the linear response approximation 
is able to describe the fidelity decay for quite a long time 
very well. Immediately beyond the time when the linear 
response formulae fail, the asymptotic results l|ll|) de- 
scribe fidelity decay quite well. We also compare to nu- 



merical random matrix simulations, where we computed 
averages over 10 4 samples of 100 x 100 matrices. Only 
the 10 central states have been taken into account. In the 
GOE case, we have also performed numerical simulations 
for symmetric perturbations with deleted diagonal (case 
(ii)) and we have not found significant deviations from 
the antisymmetric case (i) . For strong perturbations the 
plateau disappears and we find a partial revival of fidelity 
near the Heisenberg time, similar as in Ref. |l6| . 

We have concentrated on the unitary antisymmetric 
perturbation, because of its invariance properties, which 
allow to obtain results in closed form. Yet the linear 
response result in Eq. @ can be carried to higher order, 
and at least up to sixth order they coincide with the one 
for symmetric perturbations with deleted diagonal (ljj . 

The RMT model can also be compared to dynamical 
systems with chaotic classical limit. For this purpose we 
have considered a quantized kicked top [2(| ■ 

In the first example, we choose a one step propagator 



i\S x 



7 = tt/2.4 (13) 



with P = e~ laS z/ 2S e~ lSz and Sx,y,z being standard spin 
operators. Uq is time-reversal invariant, and the pertur- 
bation S y is antisymmetric in the eigenbasis of Uo- The 
"symmetrization" of Uo is essential for V to anticommutc 
with the time-reversal symmetry. 

We choose the spin S = 200, one initial random state 
and average the fidelity over 400 realizations of the prop- 
agator U\ where for each realization we draw a param- 
eter a from a Gaussian distribution of width 1 centered 
around 30. The results of fidelity decay for different 
strengths of perturbation are shown in Fig.^a). We find 
good agreement with the square of the theoretical result 
1)12(1 for the fidelity amplitude, which in turn agrees well 
with RMT simulations for the fidelity (not shown). With- 
out averaging over an ensemble of dynamical systems we 
get considerable fluctuations around RMT curves. Note 
that we do not use any fit parameters. The dimension- 
less perturbation strength e in Eq. 1)12)1 is obtained as 
e = 2Na cl {S\) 2 = 4\ 2 S 3 a ch where cr c i = 0.153 is an 
integral of the classical correlation function calculated 
using the corresponding classical map, see e.g. Ref. 
for more details. Heisenberg time is in = N = 25. 

We also consider an unperturbed propagator without 
time-reversal symmetry, that corresponds to GUE case, 



U X = Pc- i ^e- i ^/2S e -i£S Xe 



-iAS' 



(14) 



with 7 = 7r/2.4, fi — 10, £ = 1. Here we have set diagonal 
matrix elements of the perturbation in the eigenbasis of 
U to zero by hand, S' x = S x - diagS" x . We take S = 200 
and average the fidelity over 100 samples, similarly as for 
GOE case. As above we determine e from the classical 
correlation integral cr c \ = 0.168. In Fig. \Bp the results 
of the numerical simulation are plotted, together with 
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FIG. 2: Fidelity freeze for a quantized kicked top. In (a) 
perturbation is imaginary antisymmetric with unperturbed 
dynamics having antiunitary symmetry 113H while in (b) the 
unperturbed evolution as well as the perturbation have no 
symmetries left (I14H . The dashed lines give the numerical 
simulations, while the solid lines give the square of the theo- 
retical prediction 11 1 2t in (a) and RMT simulations in (b). 



RMT Monte Carlo simulation (full line). Again, good 
agreement with the RMT model is observed. 

We have presented RMT models that display the em- 
inent features of quantum freeze of fidelity under a wide 
range of circumstances not previously considered. We 
allow for any unperturbed Hamiltonian or ensemble of 
Hamiltonians for which the spectral form factor is known. 
The perturbations are represented by ensembles of ran- 
dom Hermitean matrices with zero entries on the diag- 
onal. We give a perturbative solution for the general 
model, and we present an exact solution obtained by su- 
persymmetric techniques, for the case of Hermitean anti- 
symmetric perturbations of GOE Hamiltonians. Kicked 
top models display excellent agreement with the random 
matrix results. 

The physical importance of such systems becomes ap- 
parent in two quite different aspects. On one hand mean 
field theories in some sense include the diagonal part of 



the perturbation in the unperturbed Hamiltonian, and 
thus the quantum freeze sheds new light on the surpris- 
ing success of such theories. On the other hand this result 
shows that for a quantum information process to be ef- 
fective beyond the Heisenberg time, one has to suppress 
the diagonal part of any static perturbation. 

We thank H. Schomerus for drawing our attention 
to Ref. 0|. THS acknowledges support by the grants 
No. 10803 (DGAPA) and No. 41000-F (CONACyT). MZ 
thanks the Alexander von Humboldt Foundation. TP 
acknowledges support by the grant P 1-0044 of the Slove- 
nian Research Agency. 



[1] A. Peres, Phys. Rev. A 30, 1610 (1984); R. A. Jalabert 
and H. M. Pastawski, Phys. Rev. Lett. 86, 2490 (2001); 
Ph. Jacquod et ai, Phys. Rev. E 64, 055203 (2001); 
N. R. Cerruti and S. Tomsovic, Phys. Rev. Lett. 88, 
054103 (2002); F. M. Cucchietti et ai, Phys. Rev. E 65, 
046209 (2002); G. Benenti and G. Casati, Phys. Rev. E 
65, 066205 (2002); J. Emerson et ai, Phys. Rev. Lett. 
89, 284102 (2002). 

T. Prosen, Phys. Rev. E 65, 036208 (2002). 
T. Prosen and M. Znidaric, J. Phys. A 35, 1455 (2002). 
T. Prosen and T. H. Seligman and M. Znidaric, Prog. 
Theor. Phys. Supp. 150, 200 (2003). 
M. A. Nielsen and I. L. Chuang, quantum computation 
and quantum information, (Cambridge University Press, 
Cambridge, 2000). 

T. Prosen and M. Znidaric, J. Phys. A 34, L681 (2001) 
T. Gorin, T. Prosen and T.H. Seligman, New J. Physics 
6, 19 (2004) 

R. Schafer, T. Gorin, T. H. Seligman, and H.-J. 
Stockmann, New J. Phys. 7, 152 (2005); R. Schafer, H- 
J. Stockmann, T. Gorin, and T. H. Seligman, Phys. Rev. 
Lett. 95, 184102 (2005). 

T. Gorin, T. H. Seligman, and R. L. Weaver, Phys. Rev. 
E in press, E-print: nlin.CD/0508008 (2005). 
F. Haug, M. Bienert, W. P. Schleich, T. H. Seligman, and 
M. G. Raizen, Phys. Rev. A 71, 043803 (2005). 
T. Prosen and M. Znidaric, New J. Physics 5, 109 (2003). 
T. Prosen and M. Znidaric, Phys. Rev. Lett. 94, 044101 
(2005). 

T. Gorin, T. Prosen, T. H. Seligman and M Znidaric, to 
be submitted to Phys. Rep. (2006). 

L. Viola, E. Knill and S. Lloyd, Phys. Rev. Lett. 82, 2417 
(1999). 

M. L. Mehta, Random Matrices, 2nd Edition (Academic 
Press, New York, 1991). 

H.-J. Stockmann and R. Schafer, New J. Phys. 6, 199 
(2004); Phys. Rev. Lett. 94, 244101 (2005). 
J. J. M. Verbaarschot and H. A. Weidenmiiller and M. 
R. Zirnbauer, Phys. Rep. 129, 367 (1985). 
H.-J. Stockmann and H. Kohler, in preparation. 
Y. V. Fyodorov and H.-J. Sommers, Z. Phys. B 99, 123 
(1995). 

F. Haake, M. Kus and R. Scharf, Z. Phys. B 65, 381 
(1987). 



[6] 
[7] 

[8 



[9 
[10 

[11 
[12 

[is; 
[14; 
[is; 
[ie; 

[17 

[18 
[19 

[20 



